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Introduction 



We consider pointwise semiclassical spectral asymptotics i.e. asymptotics of 
^ e(x,x, 0) as h — )■ +0 where e(x,y, r) is the Schwartz kernel of the spectral 

projector and consider two cases when schort loops give contribution above 
0(/) 1 - d ): 

(i) Schrodinger operator in dimensions 1,2 as potential V = ==>- 
(ii) Operators near boundaries. 



Chapter 1. Schrodinger operator 2 

This article is a rather small part of the huge project to write a book and 
is just part of subsection 5.2.1, appendix 5.A.1 and section 8.1 of V. Ivrii 
[Ivr2] consisting entirely of newly researched results. 

1 Schrodinger operator 
1.1 Main assumptions 

Consider Schrodinger operator 

(1.1) A = Y^Pjg jk Pk+V, Pj = hDj-Vj 

j.k 

assuming that 

(1.2) A is self-adjoint operator in i? 2 (X,H), M = C D , D = 1 for the 
Schrodinger operator, 

(1.3) B(0,1)CX, 

(1.4)^ \D a g Jk \<c, \D a Vj\<c, \D a V\ < c Va:\a\<K, 

(1.5) e <J2g Jk ( x Mk\t\~ 2 < c Vx 6 6(0,1) V£GM d \0. 

j.k 

Remark 1.1. A (unitary) gauge transformation u — > e lh ^^u with a real- 
valued function cp can be applied which for the Schrodinger and Dirac 
operators is equivalent to the transformation Vj \-¥ Vj + djip. Hence one can 
weaken condition (1.4) 2 , replacing it with 

(1.4)2 \ Da Fjk\ <c Va:\a\<K 

where 

(1.6) F Jk = d k Vj - djV k 

are components of the tensor magnetic intensity. 
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Here the principal symbol of A equals 

(1.7) a(x,0 = J>*& - Y/)(& - V.) + V 

and the subprincipal symbol vanishes and hence the coefficients in the 
spectral asymptotics are 

(1.8) K (x, r) = (2n)- d zu d (r - V) d + /2 y/g, 

(1.9) «i(x,r) = 

where ri = — oo, r2 = r, cc^ is the volume of the unit ball in IR d and 
g = det(g J ) -1 . We prefer to use notation k„(x, r) rather than k x „(t). 

1.2 General theory 

Let us consider spectral asymptotics without spatial mollification. First of 
all, let us check £-microhyperbolicity condition; it is 

(1.10) |V|>e Vxe 6(0,1) 
and we immediately arrive to 

Theorem 1.2. Let conditions (1.1), (1.2), (1.3), (1.4), (1.5) and (1.10) be 
fulfilled. Then 

(i.ii) IO < Co/7 1 -^) vxee(o^). 

Now let us use the same rescaling scheme as in the alternative proof of 
theorem 5.2.2. Then we should replace (5.2.16)* by 

(1.10)* \V\>e p 2 VxG6(0,7) 

and therefore we need to pick up the scaling functions 

112 1 

(1.12) 7 = e|V|+7, p = 72 with 7 = -h* and p = h^; 

we select the 7 equal /?3 to keep p7 > /?. Obviously |V X 7| < §• 
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We immediately arrive to the estimate 

(1.13) |RW it | < Ch^p d ~^-H0 x 

Ch^V^-W^lVl-i) as \V\>hl, 

2 h 2 / 3 2 

C/)-3 rf ^( ) as \V\< h*. 

So, we proved 

Theorem 1.3. Let conditions (1.1), (1.2), (1.3), (1.4) and (1.5) be fulfilled. 
Then estimate (1.13) holds. 

Corollary 1.4. Let conditions of theorem 1.3 be fulfilled. Let us assume 
that 

(1.14) {)(t)t~ s is a monotone decreasing function with some s > 0. 
Then as d + 3 > s estimate (1.11) holds and as d + 3 < s estimate 

(1.15) IO < Co/7K- d+s M£) VxgS(O^) 

The reader can reformulate theorem 5.1.14 for the Schrodinger operator 
and improved asymptotics with the "no loop" condition. Finally the results 
of section 4. 1 yield 

Theorem 1.5. Let conditions (1.1), (1-2), (1.3), (1.4) and (1.5) be fulfilled. 
Then for r < V* = infe(o.i) V the following estimate holds 

(1.16) |e(x,y,r)| < C'/T d (l + 7 ) -/ Vx,y G 6(0,l-e) 
where C = C'(d, c, I, e) and /, e > are arbitrary. 

Here we use that IWI < CIV— VAi. 
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1.3 Asymptotics without spatial 
mollification and short loops 

In this subsubsection we consider very special case d = 1 and d = 1, 2 when 
estimate |R)f | < Oj 1 ^ fails replaced by |R)f| < C/i~i and ft x < Ch~t as 
c7 = 1,2. More precisely, as |V| > hi estimates \R™\\ < Ch 1 ~ d \ U|( d - 3 V 2 
holds. 

Our purpose is to improve these estimates, possibly adding a correction 
term associated with the short loop. As |V X V| is our foe rather than our 
friend here, we assume that condition 

(1.17) |V| + |W| >e in 6(0,1) 

holds. Later we will get rid off it by scaling. 
Our goal is to prove 

Theorem 1.6. Let conditions (1.1), (1.2), (1.3), (1.4), (1.5) and (1.17) be 
fulfilled. Then 

(1.18) |e(x,x, 0) - K (x)h- d - h-l d \VV{x)\l d Q(W(x)h-l)g\ < 

CV^x)^- 2 ) 

where in the correction term W(x) x V(x) and Q will be defined by (1.45), 
and (1.63) respectively and 

(1.19) Q(X) = 0{\-^ d+3) ) as A^oo 

(see (1.64) for its asymptotics as A — > +00J and j(x) x |l/l/(x)| + hi, 

(1.20) \VV\ g =Q2g jk djV-d k V)K 

j.k 

The crucial step in the proof is 

Proposition 1.7. Let u(x,y, t) be the Schwartz kernel of e lh A . In frames 
of theorem 1.6 

(1.21) \F t ^ h -i T x T (t)r x u\ < Ch 1 ~ d >y(x)zV-Q V r : |r| < e 7 (x) 

where T is the small constant. 
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Note that the rescaling arguments yield (1.21) with T = eo7(x)5 and 
this would yield estimates 0(h 1 ~ d 'j( d ~ 3 ^ 2 ). Generalizing it to T = eo leads 
to remainder estimates 0(h 1 ~ d '~f( d ~ 2 ^ 2 ) which are announced in theorem 1.6, 
but with Tauberian main part. This Tauberian expression would have two 
contributors: one equal h~ d hq(x) which is the main term from t = and 
another equal to the correction term which from the loop and also collects 
all other terms from t = 0. 

Proof of proposition 1.7. Consider first the left-hand expression of (1.21) 
with Xr(t) replaced with Xr(t) supported in ^T < \t\ < T with T = Tq 
which is a small constant and with 7(x) < ei 7"; then it would be less than 
Ch s . 

Rescaling x h-» x new = xT~ 2 , t \— > t new = tT^ 1 , h \-> h new = hT~ 3 and 
multiplying operator by T~ 2 we reduce the general case of T > 7(x)z to the 
previous one; so now the left-hand expression of (1.21) with Xr{t) instead 
of Xr{t) does not exceed 

where 7" 1_2d = T x T~ 2d and the first factor T appears from Fourier 
transform while T~ 2d appears because we have a density. 

Summation with respect to T G [C07 5 , To] results in the same expression 
as / = 72 1. e. 

c-n-'iAy 

7 2 

which does not exceed the right-hand expression of (1.21). 

Therefore we need to prove (1.21) with T = C072. Rescaling x \-t x new = 
X7 _1 , t 1— y t new = t7 _ 2, h \-¥ h new = hj~2 and multiplying operator by 7 _1 
we reduce (1.21) to the case 7x1. 

2 

However, if originally 7(x) > Co/)3 then condition (1.10) is fulfilled after 

rescaling this estimate follows from (4.1.82) of [Ivr2]. On the other hand, if 

2 
originally 7 x hi then after rescaling h new x 1 and this estimate holds as 

well. □ 

Therefore due to Tauberian theorem we arrive to 
Corollary 1.8. In frames of theorem 1.6 



;i.22) \r\ vX \ < ch^^xy^-^ij) 
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as Nj L is defined with T = Tq which is a small constant. 

Now our goal is to calculate Nj with the indicated error. Due to the 
same arguments as in the proof of proposition 1.7 

(1.23) Estimate (1.22) remains true as Nj . is defined with T = max(C72, hi~ s ). 

Now let us calculate. We will do first the general calculations in the case 
of 7(x) x 1 and then we rescale. To do this we need to prove 

Theorem 1.9. Let x, y be fixed points and let £i < t 2 be of the same sign 
and t\ x t 2 x (£2 — £1) x 1- Assume that 

(1.24) There exist only one Hamiltonian trajectory (x(t),£(t)) such that 
x(0) = y, x(t) = x and t G [ti, £2]/ let it happen as t = t, £(0) = fj and 
(t- tl )x(t 2 -t)xl ; 



dx 

: l25 » * 



_ dx 
t=t ~ dt 



xl 

t—t 



and 

(1.26) Map Zo l~l TyX 3 r\ — > 7i x ty t (y,n) is nondegenerate in (y,fj). 
Then as t\ < t — e, t2>t + e,x = x,y = y 

(1.27) h- 1 ! (F t ^- lr x £ (t-7-*) ( y)c/r = e / ''^^)^b n (x,y)/) 1 

J -°° n>0 

where 

(1.28) 0(x,y)= !\{x{t),i{t))dt, £(x,0-=(d ( a,0- 

Jo 

Proof. First, let us rewrite the left-hand expression of (1.27) as 

(1.29) i^F^-i^Xeit-li 



T = 



On the other hand we know that 



;i.30) u{x, y, t) = {2nh)-^ d+m ^ f e ''^V(x,y,e,t) J^ b' n {x, y, t, 9) d6 

J n>0 
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where is m-dimensional variable and ip(x,y, t,$) is denned in the corre- 
sponding way. For example, one can take m = d, 9 = rj, 

(1.31) <p(x, y, t, rj) = — (y, rj) + ip(x, t, rj) 
where 

(1.32) d t ip = -a{x,d x ip), ip{x, 0,r}) = {x,rj). 

Let us apply stationary phase method with respect to 6, t; condition (1.24) 
and d t ip = —a(x, £) imply that there is only one stationary point and it is 
(fi, t). Further, conditions (1.25) and (1.26) imply that this is non-degenerate 
point. Thus we gain a factor hz(d + 1) and (1.27)-(1.28) are proven. □ 

Proof of theorem 1.6 To apply theorem 1.9 to our case we need just rescale 
x i— ¥ X7 _1 , h \— >■ h — /?7~5; then as K < h s i. e. 

(1.33) 7 > /)§ (1 - 5) 
we conclude that 

(1.34) N x , corr := Nj - N™ = e 1 ^ £ b n (x, x^-^-^h^^ 

n>0 

where we also multiplied by 7~ d since we are dealing with densities. 
Thus 

(1.35) Correction term N xcorr is of magnitude h^~ d '^~^ d+ ^> . 

This implies a drastic difference between d = 1 when the correction term is 
below remainder estimate h 1 ~ d/ y(x)2^ d ~ 2 ^ of (1.22) only as 7 x 1 and d > 2 
it is always so as 

2d-2 

(1.36) 7 > 71 := /i33=i 

Problem 1.10. Consider averaged with respect to spectral parameter correc- 
tion term and to prove that it is of magnitude 

/j |(l-d) 7 -J(«/+3) 7? /A.J < ^(1-^-1(^3+25^ A 

L72 L 

under condition 
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(1.37) 1l)(t)t~ s is a monotone increasing function with some s > 0, 

I believe that the complete proof of this statement is worth to be pub- 
lished. 

Let us introduce 

(1.38) X ={x: V(x) < 0}, X° = dX~ = {x : V{x) = 0}, 
Without any loss of the generality one can assume that X~ = {x 1 < 0}. 

1.3.1 Case d = 1. 

We can assume without any loss of the generality that a(x, £) = (3(x) (& — 
Vo(x)) (we can always get rid off V\ by gradient transform); then 

Xl 



;i.39) l^(x 1 )=(^| Xl \/o(y 1 )^y 1 



is a travel time from xi G X to X° (on energy level 0, if we replace /3 by 1) 
and one can see easily that then 

(1.40) Under assumption (1.17) l/l/(x)/Vo(x) is a smooth and disjoint from 
function on X U X°. 

We redefine x\ = 1/1/ (x) and then a(x, £) will be in the same form as before 
but with different /3 and with l/l/(xi) = xi: 

(1.41) a(x,0 = /3W(e i 2 -x 1 ) 

Let us prove that in calculations one can replace (3(x) by (3(x) = 1 and 
assume that 

(1.42) A = h 2 D\-x y on R. 

Really, let x be a point where calculations are done, while x be a "running" 
point. Without any loss of the generality one can assume that 0(x) = 1. 
Let us rescale as before. Then we can assume that we are at he point with 
7(x) but on [— c, c] 

(1.43) \D{f3(x)\ < qel 
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where e is an original 7(x). 

Let us apply theorem 1.9. Note that phase functions for the original 
operator and for the model operator coincide identically while amplitudes 
differ by 0(e) (where e = 7); so an error is 0(h 1 ~ d e) and scaling back we get 
an error estimate Ch 1 ~ d ^2( d ~ 3 )e = Ch 1 ~ d ^2( d ~ 1 ) which is actually better by 
factor 72 than we need. Calculations for model operator (1.42) are produced 
in Appendix 1.4. Theorem 1.6 is proven as d = 1. 



1.3.2 Case d > 2. 

Again without any loss of the generally one can assume that V(x) 
with k(x) > disjoint from 0. 
Let 



-k(x)x t 



(1.44) = {(x, : V(x) = 0}nl„ = {(x, £) : x x = 0, £ = Vj(x)} 

parametrized by x' = (X2, ... ,x</). 

Consider Hamiltonian trajectory passing through (0,x') G 0. One 
can select £1 as a natural parameter along this trajectory, so we get a d- 
dimensional manifold A C 5Z ; since Zq is (2d — l)-dimensional, there exists 
(d — l)-dimensional variable 77: A = Zo fl {rj = 0}. 

Assume first that 



(1.45) There is no magnetic field, i. e. V 1 = V 2 



V H = 0. 



Then trajectory passing through as t = is symmetric: x(— t) = x(t) and 
£(— t) = (,(t). Then the set of loops coincides with A and x-projection of 
the loop is exactly as on the left picture below 



x(0) 



x(t)=x(-t) 



*(0) 



■x(t) = x(-t) 



(a) Original coordinates (b) Straighten coordinates 

Figure 1: : Short loops 
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If we introduce new coordinates replacing x' by x'(0) we get picture as 
on the right. Note that due to the same arguments as for d = 1 one can 

assume that 

(1.46) a(x, = /*(*) (£ - *i + X^^Wfe-^MeOfe-a/tMei)) 

j,k>2 

with positive definite matrix (g' J ). Then 77 = if and only if f — a/(x)£i = 
V/ = 2, ... , d. Note that on Zq we must have i] = ==>- {a, r^} = i. c. 
{£i — xi, £y — aj(x)^i} = which easily yields that a: 7 (x) — V/ = 2, ... , d i. 

c. 

(1.47) a(x, = /?(x) (e? - X! + £ ***(*)&&) 



j,k>2 



and r/ = f . 



x(0) 



One can see easily that in the special case 
g'-' k (x) = 8jk trajectories are parabolas 




6 
f 



const, 



Xl 



X 



const + ft (y > 2) 



Figure 2: Selecting f ^ de- 
stroys the loop 



and similarly looking in more general 
case. 

Again, as 7 x 1 we can estimate con- 
tribution of the loop by Ch 1 ~ d x h^ d ^ = 
/j^z^ -1 ) where extra factor h^ d ~ x > due 
to theorem 1.9. Scaling procedure leads 
to estimate (1.35). 
As g'J k = Sjk Appendix 1.4 implies that estimate (1.18) holds. Thus it 
holds for g'J k = const. For variable g'J k we get instead of phase |f | another 
phase |f | +ou(x, f ) with ui = 0(7 |f | ) which leads to an error with an 
extra factor 7 in the error estimate which becomes /7 _ 2( ' _1 )^~4( d_1 ) and is 
less than the right hand expression of (1.18). 

So, under assumption (1.45) theorem 1.6 is proven. 
We can consider operator (1.47) and add magnetic field to it. 

In this case equalities £(— t) = — £(t) 
and x(— t) = x(t) along trajectories pass- 
ing through fail, A loses is value and 



<(0) 



Figure 3: Short loop with mag- 



nci"i _ ir' fiolri 
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picture on figure 1(b) is replaced by 
the picture on the left where actually 
the arcs are symmetric only asymptoti- 
cally. However as before for each x (with 
< xi < e) there exists £ = £(x) such 
that (x, £) G Zo and Hamiltonian trajec- 
tory passing through (x, £) comes back to 
x. So, instead of few looping trajectories 
and each such trajectory being a loop for 

each of its points, we have now many looping trajectories but each of them 

serves only one point. 

Further, one can prove that then the phase is changed by 0(e) where 

3 i 5 

e = 7 and scaling back we get an extra term 0(h~ 72 x e) = 0(h~ 72) 
in the exponent which leads to the error not exceeding Ch~^ d+1 '^~^ d ^ 7 ' 

2 

which does not exceed the right-hand expression of (1.18) as hi < 7 < 71. 
Theorem 1.6 proven completely. □ 

Remark 1.11. Obviously 

4 3 2 

(1.48) W(x) = const Q*{x), const=(-) 3 

where q(x) is the distance in the metrics g-' k V(x)^ 1 from x to X°. Thus one 
can rewrite the correction term as 



(1.49) /ri d |V\/(x)|5g(const/rM( x ) 

Finally, let us get rid off condition |W(x)| > e. To do this let us rewrite 
first the right-hand expression of (1.18) as 

(1.50) C/r3( d - 1 )(max(|\/|,/?§))^ 2) 

thus releasing notation 7. 

Let us introduce a standard scale 7 = (| V| + 1 V V| 2 ) 2 + h? but exclusively 
at point x in question. Then applying the standard rescaling x — y X7 _1 , 
h -»■ h = h-f~ 2 , V -»■ Vj- 2 , W -> W7- 1 we get instead of (1.50) 
expression 

C/) 1 -V" 2 (max(|\/|7- 2 ,/ J i7-t))^ d - 2) = CM- d >(max(| V|, /iM))^ -2 *. 



'1.52) C< 
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As d > 2 the latter expression can only increase as 7 is replaced by 1. 
However as d — 1 situation is different and we arrive to 

(1.51) CmindUI^/r^-s)). 

Thus we proved 

Theorem 1.12. Let conditions (1.1), (1.2), (1.3), (1.4) and (1.5) be fulfilled. 
Then 

(i) As d = 2 asymptotics (1.18) holds; furthermore one can skip the correction 
term without penalty as h d ~ 1 \W\ d+1 < IVp^ -1 ; 

(ii) As d = 1 the left-hand expression of (1.18) does not exceed 

' \V(x)\~s as /j§|W(x)|i < \V(x)\ < |W(x)| 2 

or as \V(x)\>max(\W(x)\ 2 ,h), 

/T3|W(x)|-3 as \V(x)\ < /?i|W(x)|i, \VV(x)\>h*, 

h- 1 * as \V(x)\ < h,\VV(x)\ < /ii; 

further, as |V(x)| > max(|W(x)| 2 , /?) one can skip the correction term 
without penalty. 

1.4 Spectral kernel calculations for some 
model operators 

1.4.1 Case d = 1 

Consider operator (1.42). 
Note first that 

(1.53) B = -xi ==>• 9 r ee(x,y,r) = 6(xi -yi)6(xi +r). 
Then making unitary /?-Fourier transform we conclude that for 

(1.54) B = /i% =► 9 r e s (^i,77i,T) = ^^exp^"^ - ^)r). 
Therefore since T*hD^T = hD^ + £ 2 for T = exp(//? -1 |£i) we conclude that 

(1.55) e = /!% + £ =» 

9 T e B (6, »7i, r) = (27T/))- 1 exp (/fT 1 [(& - ^)r - ^ 3 ~ Vi)] ) ■ 
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Finally making inverse Fourier transform we conclude that for A denned by 

(1.42) 

(1.56) <9 T e(xi,yi,r) = 

(2nh)- 2 J J exp(//)- 1 [(& - m )r + x^i - ym -\(ti- vl)]) ^i cfyi 
and therefore 

(1.57) <9 T e(xi,xi, t) = 

(2nh)- 2 J! expiih- 1 ^ - m )(r + Xl ) - I^l-Vl)}) d^d m . 
One can rewrite it as 

_ 2 

/(xi+r)h 3 
F{r)dr, 
-oo 

with 

(1.59) F(t) = (2vr)- 2 ||exp(/[(6 - 77i)t - ^ 3 - ^ 3 )]) ^i d m = 

2(2tt)- 2 f f exp(i[2fit - ^fi 3 - 2fia 2 ]) da dp = 

(27r)-i y |^|-i exp(-/^ sign /3 + i[20t - ^/3 3 ]) c/£ 

where we first substituted £i = a + (3, r\\ = a — j3 and then integrated with 
respect to a. 

The remaining integral has stationary points /3 = ±t2 and a singular 
point /3 = 0. Decomposing for t > F(t) = Fi(t) + F 2 (t), 

(1.60) Fi(t) = (27r)-i / |/3|-5 exp(-/^ sign /3 + 2//3t) t//3 = ^«H 
with 

/■oo 

(1.61) /t = 2(27r)-i / /H cos(2/3 - -) c//3, 

Jo 4 

and 

(1.62) F 2 (t) = 

(27r)-i 1 1/5| i exp(-/^ sign /3 + 2//3t) (exp(-^/3 3 ) - l) c/0. 
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Then primitives of Fi(t) and /"2(f) produce exactly main term and 
correction term in (1.18): 

(1.63) Q{t)= f F 2 {t')dt' = 

- 2- 2 (2vr)-i J \0\-l (/£ sign /3 + 2i/3t) (exp(-^/3 3 ) - l) d/3. 
In this integral singularity at /3 = gives a relatively small contribution 

3 

(one can prove it is 0(t~z) as t — > 00); the main contribution comes from 

1 3 

the stationary points (3 = ±£2 and modulo 0(t~?) 

(1.64) Q(t) = -2- 1 (2tt)- 1 £ |/5| 2 exp(/^sign/5 + ^//5t) = 



p=±ti 



(27r)- 1 t- 1 sin(^i). 



1.4.2 Case d = 2 

Now we consider operator 



A = h 2 D 2 Xi + h 2 D 2 X2 - Xl 



leading to B = £ 2 + £ 2 + /?D^ and then instead of (1.57) we have 
(1.65) 9 T e(x 1 ,x 1 , r) = 

(2nh)- 3 J J exp(//)- 1 [(& - ^)(r + x, - £ 2 2 ) ~ ^l - v!)] ) d^d^dm- 

and all previous formulae are adjusting accordingly. Then (1.58)-(1.59) are 
replaced by 

_ 2 

/(xiH-r)ft 3 
F{r)dT, 
-00 
with 

(1.67) F(t) = 2(2tt)- 3 /"/"exp(/[2/3(t + 7 2 )-^/3 3 -2/3a 2 ])c/«c//3c/7 = 
(27r)- 2 y |/3|- 1 exp(/[2/3t-|/3 3 ]) C //3 
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and (1.60)- (1.62) are replaced by 

(1.68) F^t) = (2tt)- 2 / \fi\- 1 exp(2//?t) dfi = n + Kl r 3 
with 

/■OO 

(1.69) K = 2(27r)-t / p-zcos{2j3--)dj3, 

Jo 4 

where integral was understood in the sense of distributions and 

(1.70) F 2 (t) = (2n)- 2 r 1 J |/3|- 1 exp(2//3t)(exp(-^//3 3 ) - l) c//3. 
Then (1.63) is replaced by 

(1.71) Q(t)= f F 2 {t')dt' = 

J OO 

-2- 1 (27r)- 2 y'|/3|- 2 sign/3exp(2//3t)(exp(-^//3 3 )-l)d/3. 

In this integral singularity at /3 = gives a small contribution; really, 
decomposing exp(— 1//3 3 ) — 1 into powers of (3 3 we get that he leading term 

is 0(r~ 2 ). The main contribution comes from the stationary points (5 = ±£2 
and modulo 0(t~2) 

(1.72) Q(t) = (27r)- 1 rfsin(^tl). 

2 Spectral Asymptotics near Boundary 

In this section we consider boundary layer type term in e(x, x, r) appearing 
near boundary. Let us recall that due to rescaling technique 

(2.1) e(x,x,r) = N^ + 0(/) 1 -V 1 W) 

where 7(x) = dist(x, dX) and this remainder estimate generated after inte- 
gration 0{h 1 ~ d \ log h\) and basically the whole chapter 7 of [Ivr2] we spent 
to eliminate logarithmic factor. 

Now we would like to consider e(x, x, r) without integration with respect 
to x. 
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2.1 Preliminary analysis 

2.1.1 Discussion 

Let us recall that in the previous chapter we derived asymptotics of 

(2.2) N:=r((? lx e(.,.,r) t (? 2y ) 
and 

(2.3) N / :=r / g x g y ((? lx e(.,.,r) t (? 2y ) 

with the remainder 0(h 1 ~ d ) provided problem (A, 56) is microhyperbolic in 
multidirection T = (£', v\, ... , vm)\ further, we derived asymptotics of 

(2.4) N x :=r x 9 x g y ((? lx e(.,.,T) t (? 2y ) 
and 

(2.5) N x :=r x ((? lx e(.,.,r) t (? 2y ) 

with the remainder 0(h 1 ~ d ) provided £' = (^,0) and for (2.5) we also 
need to assume that ui = ■ ■ ■ = i>m in the microhyperbolicity condition; 
this asymptotics had a boundary layer type term h 1 ~ d T(x',xih~ 1 ) with 
D a T(x', r) = 0{r-°°) as r ->• +oo. 

Let us discuss this microhyperbolicity condition. Consider first the case 
of the Laplace operator; without any loss of the generality one can assume 
that a(x, £) = a'(x, £') + ^\ where a(x, £') is a positive definite quadratic 
form. Then microhyperbolicity means that 

(2.6) (4,V^') + 2^i>0 as £i = ±(r - a(x, £'))* 

and we can always take l't = £', v+ = v_ = as long as a'(x, £') > (i.e. 

evo). 

So, with the exception of the case £' = both microhyperbolicity 
conditions hold. In this exceptional case however VI a' (x, £') = and the 
first term in (2.6) is no matter what £^ we pick up; so (2.6) boils up to 

±v ± (T-a(x,?))*>0 

and for microhyperbolicity we just take v± = ±1 provided r > 0; this is the 
standard ^'-microhyperbolicity condition (we deliberately take £' = (£1,0)). 
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However it is not the case as we assume that v + = v-\ then condition is 
impossible to satisfy. So, the previous section fails to find asymptotics of 
e(x, x, r) but not because of the rays tangent to the boundary but on the 
contrary, because of the rays orthogonal to it (in the corresponding metrics). 

This is a pleasant surprise because then we can find the solution of 
non-stationary problem in the form of the standard oscillatory integrals 
as we did in section 1.3. The analogy does not stop here: the exceptional 
rays hit the boundary orthogonally, reflect and follow the same path as 
before forming the short loops of the length 2dist(x, dX) where distance is 
measured in the corresponding Riemannian metrics. 

However there is a difference: in the former case loop appeared because 
trajectory went "uphill", lost velocity and rolled back "downhill' repeating 
the path while now it reflects without losing velocity. 

2.1.2 Pilot model 

Consider X = IR + x IR d_1 , a(£) = |£| 2 and either Dirichlet or Neumann 
boundary conditions. Then propagator is defined by 

(2.7) U(x, y, t) = U°(x, y, t) + U\x, y, t) 
with free space solution 

(2.8) U°{x, y, t) = {2nh)- d f e / ^ 1 (l?l 2 t+<x-y,5» d £ 
and reflected wave 



(2.9) U\x,y, t) = ^(2nh)- d f, 



e ih- 1 (\(\ 2 t+{x-~y,0) d £ 

where y = (— y\,yi, ■■■ ,yd) and <^ = =pl for Dirichlet or Neumann boundary 
conditions respectively. 

Since U 1 is responsible to the difference e 1 (x, y, r) between e(x, y, r) and 
its free space counterpart e°(x,y, r) we manipulate only with it. Taking 

y = x we get (x — x, £) = 2xi£i while F t ^ h -i T replaces e lh T by hb(r — |£| 2 ); 
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:■;«) 



(2.10) F^,-!^ 1 = ,(2nhY- d J8(r- l^e 2 * -1 ** </£ = 

, . „ , f / cos(2/)" 1 r5xicos0)sin d ^ 2 0o'0 d > 2, 

<?const ■ h 1 ~ d T2- 1 ■ < Jo 

{ cos^/rVsxi) </ = 1 

where we introduced spherical coordinates (0,sin0£) on E> d ~ 1 , <p £ (0,7r), 
CG§ d - 2 . 
It leads to 

(2.11) e^x.x.l) = h- 1 f (F^-iSxU 1 ) dr = h^T^h^x,) 

J — oo 



with 

(2.12) T ? (r) 



const • c ^ 



/»1 /»7T 

/ z^ 1 / cos(2zrcos0)sin d ^ 2 0c/(/.cyz cf > 2, 

Jo Jo 

cos(2zr) d = 1 



o 

2 



where we plugged t = z . 

Here we took r = 1 since any other value could be reduced to this one 
by rescaling h i— >■ hr~z. 

One can find the constant so that for Dirichlet/Neumann problem 
e 1 (x, x, r) = <;e (x, x, r) as Xx = 0. 

Stationary phase method implies that 

(2.13) T(r) = const •r" (d+1)/2 cos(2r) + 0(r- (d+3)/2 ) as r -> +oo. 
Therefore 

(2.14) With an error not exceeding 

(2.15) Ch- d {x 1 h- 1 + 1)-^+!) + C/) 1 ^ 
the standard Weyl expression for e(x, x, r) holds; 

this error is better than 0(h~ d (xih~ 1 + l) ) which was due to rescaling 
but it is only because critical points of a(£') = |£'| 2 are non-degenerate. 
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2.2 Schrodinger operator 

Let us consider Schrodinger operator with the symbol 

(2.16) a(x,0 = J2s Jk & ~ V/)(& - ^) + V 

assuming that r — V x 1 and positive. We want to study trajectories on 
energy level r only and we can assume without any loss of the generality 

that 

(2.17) r = 0, V(x)<-e, V(x) = -1 

where the last assumption is achieved by division a(x, £) by — V[x) which 
does not affect trajectories on energy level 0. 

2.3 General theory 

We can assume without any loss of the generality that Vi = 0. Let us 
introduce coordinates xi = dist(x, dX) in the metrics (g^ k ) 1 - ) and x' which 
are constant along trajectories of a(x, 6) on level which are orthogonal to 
the boundary, i.e. such that 

(2.18) e = 6;i. 

We are interested in the construction of parametrix in 

(2.19) U = {\Zj-Vj\<€ 7 = 2 d}. 

It is be much simpler to consider this problem with "time" Xi and with t as 
one of spatial variables. Then as / = 

(2.20) u'(x, y,t)= fj2 e ih ~ 1 ^ x ' y - v)+trio) J^ f U x < ^ ^) h "~ d dr l~ 

where 77 = (77', rjo) = (772 , - - - . Vd, Vo) and phase functions 0° are defined from 

(2.21) 9 xl 0( = A,(x,V x ^ ? ,?7o), 

(2.22) <f,{x,y,n,n Q )\ x ^ yi = {x'-y',n'), 

(2.23) A,(x,e / ,77o) = <?7o-a(x,e / )) 1 . 



^ I.e. in the final run in (g jk (— V(x)) 1 ). 
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Here amplitudes f® n {x,y,n) are smooth, satisfy transport equations and 
some initial conditions as x\ = y\ but we do not need to calculate them. 
Then we can calculate u°, hDiu as x 1 = and define u 1 (x,y, t) in the 

same way (2.20) but with 0j satisfying (2.21) and 

(2.24) <f>\ = -(fp_. as X! = 0. 

Here amplitudes f* n (x,y,n) are smooth, satisfy transport equations and 
some initial conditions as x\ = but we do not need to calculate them. 
Obviously 

(2.25) 4i(y,y,v)=yLfa(y,v) as x 1 =y i , 



(2.26) ^(y, n) = \(0, y' , rj) - A_,(0, y\ n) + 0( yi ). 

Therefore (as \t\ < T) 

(2.27) r y (Q lx u 1 {x,y,t) t Q 2y ) = 

J J2 e^y^^^ J2 C(y, v)h n ~ d dn, 

•* ?=± n>0 

where Qj = Q(x, hD', hD t ) are cut-offs with respect to £', r and finally 

(2.28) F t ^ h -i T x T (t)(r y {Q lx u 1 (x,y,t) t Q2y)) = 

J2 ! e ''"-Vi^(yV.-) J2 C{y, rf, r)h 1+ "- d dn' . 

q=± •* n>0 

Obviously, representation (2.28) implies that 

(2.29) For small constant T expression (2.28) does not exceed Ch 1 ~ d . 
Since (2.29) is true for u under £-microhyperbolicity condition 

(2.30) \r-V(y)\>e 

due to free space results, this statement is true also for u and therefore due 
to Tauberian arguments 
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Proposition 2.1. Let (2.30) be fulfilled and Q 1 = Qi(x, /?D'), Q 2 (x, /jD') 6e 

operators with symbols supported in U. Then T y ((5i x e 1 (x,y, r) f (?2y) equals 
modulo 0(h 1 ~ d ) Tauberian expression 

(2.31) h- 1 f T (F M ,-vXr(t)(r y (Qix« I (x,y,t) t (3 2y ))) dr'. 

J —CO 

Due to results of subsection 7.3.1 of [Ivr2] we need to calculate only 

(2.31) with an extra factor y(r — r') with y supported in the small vicinity 
of 0: 

(2.32) f X (r - r') E / e^ 1 **^^ £ £(y. »/. r')/^ <V<^'. 

Obviously modulo 0(/7 1_c/ ) we can skip in the right-hand expression of (2.32) 
all the terms with n > 1; so we get h~ d J with obviously defined J: 

(2.33) J := J2 f X(r ~ r') f e^^'^f^y, rf, r') dn'dr' 

with h = hy^ 1 . 

Further, for y\ > h^ 1 one can consider J as an oscillatory integral 
with a semiclassical parameter h~. One can see easily that due to (2.25)- 
(2.26) in this oscillatory integral the phase function Vv(y> ??> r ') satisfies 
l^i-'VvCy. v' < T ')\ > e i an( i has non-degenerate critical points with respect to 
rf; we denote these points by rf^y, r'). Therefore 

(2.34) Js^^e rI ^ T V ? ;(y,r)# 1+ ^ as fr«l 

where 

(2.35) i P,{y,r)=i/j,(y,r] , (y,T),T). 
Then 

(2.36) -4 = J° + 0(/?) where J° is also given by (2.33) but with y replaced 
by(0,/). 

Also the fact that /q are defined through transport equations along trajecto- 
ries implies then 

Proposition 2.2. Modulo 0(h 1 ~ d ) term (2.31) coincides with the same 
term but calculated for (A, 36) = (a(0,y', hD x ), 56(0,/, hD x )). 
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2.3.1 Calculations and main theorem 

Now in virtue of propositions 2.1, 2.2 we need only to recalculate T(r) in 
the case of more general boundary conditions than Dirichlet and Neumann, 
and then make change of variables. 

Repeating arguments of subsection 2.1.2 one can prove easily that 

Proposition 2.3. Consider Schrodinger operator with g jk = bjk, V = —1, 

r = and the boundary condition 

(2.37) (e (x, hD')hD 1 + iB t {x, hD')) u| xi=0 = 
with 

(2.38) \b {x',e) + \b 1 (x',e)\>e. 

Assume that supp qj C U defined by (2.19). Then 
(i) The boundary layer term is h~ d T(h~ 1 xi) with 

(2.39) T(r) = 

(2vr)- d /( f (? ' T) (£i + im')) (6 - imT 1 ^ 1 ^iW)*(0 <% 

where bj are principal symbols of Bj and due to self-adjointness assumption 
f3 = —ibi : bo is real-valued. 

(ii) Here function T(r) admits decomposition 

(2.40) T(r) ~J2e 2ir "J2 x <" rUd+1) ~ n as r ^ +°°- 

S=± n>0 

Therefore changing variables and using asymptotics in the zone where 
theorems 7.3.10 and 7.3.11 of [Ivr2] work we arrive to 

Theorem 2.4. Let us consider Schrodinger operator with self-adjoint bound- 
ary condition (2.37) and let condition V(x) < —e, be fulfilled. Further, let 
ellipticity condition 

(2.41) \b (x',£')\(x', £') + M*'.OI/0 on supp qriLI supp q 2 
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be fulfilled where X(x' , £') is a root of a(x, £i, £') = with Im A > 2 \ 

Then 



(2.42) r x ((? lx e(x,y,0) t (? 2y ) - K Q (x)h- d - h- d T(x' , h~ 1 x 1 )^g\ < Ch 



,1-d 



wi/i IVey/ coefficient k and with g = det(g j7c ) -1 where x 1 is the distance 
from x to dX in the metrics gi k \J~ x anc [ T(r) is defined by (2.38) and 
satisfies (2.40). 

Remark 2.5. (i) Due to theorems 2.4 and 2.6 the above theorem but with 
T(r) replaced by T(r) + T^r) remains true if ellipticity condition (2.41) is 
replaced by microhyperbolicity condition 

(2.43) |6 (x',OA(x',£') + M*'-£')l+ 

|V^(6o(x / ,OM* / .0 + M* / .0)l^° on suppgiUsuppq-2 
where T fa (r) G J^(R+); 

(ii) Standard Weyl asymptotics holds with the remainder estimate 0(h 1 ~ d ) 
for e(x, x, t) as xi > /,( rf - 1 )/( rf + 1 ). 

2.4 Generalizations 

Now we need to consider the generalizations to the systems. We consider 
only Qj with symbols supported in U. The same approach leads to the 
following 

Theorem 2.6. Let (A, 36) be self-adjoint and letU C T*dX. Assume that 

(2.44) As (x, £', r) G [0, e] x W x [— e, e] a// t/je roots £i o/ equation 

(2.45) det(r - a(x, £)) =0 

are A- (x,f,r) with j = 1,...,M ± which are real-valued, smooth and of 
constant multiplicities and 

(2.46) ± d T Xf > 0; 



2 ) So, this is the condition only in the complement of nll\gx where E = {(x,£) £ 
T*X,a(x,£) = 0}. 
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Then 

(i) Reflected solution has decomposition 

(2.47) F t ^ h - lT XT(t){r x (Q lx u 1 (x,y,t) t Q 2y )^ = 

J- J e'h-^^M.r) J2 f Jkn {x. rf, r)h 1+n - d dr/. 

l<J<M+,l<k<M- n>0 

with real phase functions ifijk and does not exceed Ch 1 ~ d (l + h~ 1 yi)~2^ r ~ 1 ); 

(ii) Asymptotics (2.42) holds with h~ d T(x, h~ 1 xi) constructed from (2.47) 
in the usual Tauberian way; thus 

(2.48) T(x,/?- 1 x 1 ) = 

l<J<M+,l<k<M- f)>0 

(Hi) Under assumption 

(2.49) r Vj(A+(x,£',r)-A*(x,£',r))=0 

=► rankHess^(A+(x,^,r)-A fc (x,^,r)) >r 

expressions (2.47) and (2.48) do not exceed Ch 1 ~ d (l + h~ 1 x{)~^ r ~ 1 ^ and 
C(l + /)~ 1 x 1 )~2"( r + 1 ) respectively; 

(iv) Under condition (2.49) 3 one can replace x by x' in the phase functions 
and amplitudes in (2.48); then T(y', /j _1 yi) is constructed for the problem 
(a(0,y',hD x ),db(y',hD x )); 

(v) As r = d — 1 h~ d T(h~ 1 xi) admits stationary phase decompositions 
(2.50) J2 e ' Xlh 1>1JkiX '' T) J2 f J'<n( x '- T )( x ^ 1 y 1 ' {d+iyn mod °( hl ~ d Y 

j,k n 

Remark 2. 7. (i) An easy proof is left to the reader. 

(ii) In contrast to the Schrodinger operator we do not define xi as the distance 
to the boundary; eliminating Xi from the phases ipjk without deteriorating 
estimate is possible as r > 3 while eliminating xi from amplitudes is possible 
as r > 1. 



Bibliography 26 

Bibliography 

[Ivrl] V. Ivrii. Microfocal Analysis and Precise Spectral Asymptotics, 
Springer- Verlag, SMM, 1998, xv+731. 

[Ivr2] V. Ivrii. Microfocal Analysis and Sharp Spectral Asymptotics, in 
progress: available online at 
http://www.math.toronto.edu/ivrii/futurebook.pdf 

Department of Mathematics, 
University of Toronto, 
40, St. George Str., 
Toronto, Ontario M5S 2E4 
Canada 

ivrii@math.toronto.edu 
Fax: (416)978-4107 



